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We study the evolution of the distribution of eigenvalues of a N x N matrix subject to a random 
perturbation drawn from (i) a generalized Gaussian ensemble (ii) a non-Gaussian ensemble with a 
measure variable under the change of basis. It turns out that, in the case (i), a redefiniton of the 
parameter governing the evolution leads to a Fokker-Planck equation similar to the one obtained 
when the perturbation is taken from a standard Gaussian ensemble (with invaraiant measure). This 
equivalence can therefore help us to obtain the correlations for various physically-significant cases 
modeled by generalized Gaussian ensembles by using the already known correlations for standard 
Gaussian ensembles. 

For large iV- values, our results for both cases (i) and (ii) are similar to those obtained for Wigner- 
Dyson gas as well as for the perturbation taken from a standard Gaussian ensemble. This seems 
to suggest the independence of evolution, in thermodynamic limit, from the nature of perturbation 
involved as well as the initial conditions and therefore universality of dynamics of the eigenvalues of 
complex systems. 

PACS numbers: 05.45+b, 03.65 sq, 05.40+j 



I. INTRODUCTION 



Many important physical properties of complex quantum systems can be studied by analyzing the statistical response 
of eigenvalues to an external perturbation. It turns out that the evolution of the eigenvalues of the system as a function 
of the strength of the perturbation is universal but only after an appropriate normalization of the perturbing potential 
[1,2]. Various approches have beeen adopted for a better understanding of this phenomena e.g. non- linear a model 
[3], random matrix theory (RMT) [4], semi-classical techniques [5]ctc. Among these, a special class of Gaussian 
ensembles of RMT have been used very successfully in modelling the short energy-range behaviour of eigenvalues [4] . 
The probability density of a matrix in this ensemble depends only on the functions invariant under change of basis 
e.g trace of the matrix thus making it easier to calculate the eigenvalue distribution. However it has been conjectured 
that the local statistical properties of a few eigenvalues of a random hermitian matrix are independent of the details 
of the distribution of the matrix elements apart from a few broad characterizations such as real symmetric versus 
complex hermitian etc [4] . Therefore the success of RMT in modelling the energy-level behaviour of complex systems 
should not be restricted only to the standard Gaussian ensembles. 

In this analytical work, we study this conjecture by examining the dynamics of eigenvalues of a quantum hamiltonian 
under a random perturbation belonging to a class of (i) a general Gaussian ensembles and (ii) a non-Gaussian ensemble 
with logarithm of the distribution given by a polynomial. In both cases, the probability density of a matrix is chosen 
to contain the basis-dependent functions and therefore is no longer invariant under the basis-transformation. By 
integrating over all perturbations from such an ensemble, we show that, for large N, the probability density of 
eigenvalues evolves according to a Fokker-Planck (F-P) equation similar to the one proposed by Dyson for Wigner- 
Dyson gas [6]. The strength of the perturbation acts here as a time like coordinate. For finite N, this analogy seems 
to survive only when the dynamics is considered as a function of the perturbation strength as well as the variances of 
the matrix elements of the perturbation. 

The statistical response of eigenvalues to a changing perturbation was studied by Simons and coworkers [1] by 
using a different approach, namely, non-linear a model. By considering a hamiltonian H = Hq + xV with V as the 
perturbation and x its strength, they explicitly calculate the autocorrelation function of the energy eigenvalues for 
disordered systems within the so called Gaussian approximation (or zero mode aproximation) . Other results pertain 
to the numerical study of perturbed quantum chaotic systems, such as the chaotic billiard with a varying magnetic flux 
[7]. As shown in their work, the universality manifests itself among the density correlators only after the normalization 
of eigenvalues Aj by mean level spacing and that of perturbation by rms "velocity" of eigenvalues < (dXi/dx) 2 >; by 
universality one means here the independence of the bulk correlators in a many body system from boundary effects. 
They also conjectured the asymptotic correspondence of the two-variable parametric density correlators of complex 
quantum systems with the corresponding time-dependent particle-density correlators for the Sutherland hamiltonian 
[8]. The connection between the static correlators (for a given parameter value) and the ground state correlators of 
Sutherland is already well known [1] . 
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This correspondence was analytically shown in a recent work by Narayan and Shastry [9]. They studied the 
evolution of the distribution of eigenvalues of a N x N matrix H subject to a random perturbing matrix V taken from 
a standard Gaussian (SG) ensemble and showed that it leads to the F-P equation postulated by Dyson. Within this 
model, they proved the equivalence between the space-time correlators of the ground state dynamics of an intcgrablc 
1-D interacting many body quantum system (the Calogero-Sutherland system) and the second order parametric 
correlations of H. This analogy follows because the Dyson's F-P equation is equivalent, under a Wick rotation, to the 
quantum mechanics of Calogero model which, in thermodynamic limit, has identical bulk properties as the Sutherland 
system [9,10]. 

Our present work extends this analogy further to a more general model of hamiltonians H . We proceed as follows. 
For the clarification of our ideas, we first study, in section II, the distribution of eigenvalues P(p,r) for an arbitrary 
initial condition H and the perturbation V taken from a generalized Gaussian ensemble (later referred as GG case). 
We obtain a partial differential equation governing its evolution which, after certain parametric redefinitions, turns 
out to be formally the same as the Fokker-Planck equation govening the evolution of Wigner-Dyson (WD) gas. As the 
solution of F-P equation as well as various correlatios in the latter case are already known for many initial condiitions, 
this helps us to obtain P(p, r) as well as the correlations in the GG case, as given in section III. The section IV deals 
with the evolution of eigenvalues under a more general perturbation. We are motivated to study this case because the 
claim about universality of the dynamics can only be made if the equivalence is proved, at least in thermodynamic 
limit, for a V with distribution of more general nature. We conclude in section V which is followed by an appendix 
containing the proofs of some of the results used in this paper. 



Our interest is in the evolution of eigenvalues of H under a random perturbation V with initial state H as in 
Ref.[9]. Here both H and V are hermitian and can be assumed to have same mean spacing without any loss of 
generality. To ensure the same mean spacing for H also, we need to modify the parameterization of our perturbation 
strength x [9,10,11] from H(x) = H + xV to H(x) = H Cos(Qx) + VQ^Siniflx) with fi oc -!=. Since fl -> for 

large N, the reparameterization of the perturbation strength is inconsequential for finite x. The distribution p(V) of 
matrix V is still chosen to be a Gaussian, p(V) = Cexp(— X)fc<z a klVkl) with C as the normalization constant. So far 
we are in the same situation as in Rcf.[9] but now we depart and allow the variances for diagonal and off-diagonal 
matrix elements to be arbitrary (later referred as GG case) unlike the SG case where variance of the diagonal ones 
was twice of the off-diagonal ones. 

As indicated in [9] , the evolution of eigenvalues of H as a function of x for any perturbation V can be expressed 
as a set of first order differential equations. They turn out to be analogous to the equations of motion of a classically 
integrable system with 2N + (3N(N — 1)/2 degrees of freedom (with f3 a symmetry dependent parameter). In principle, 
the distribution of eigenvalues can be obtained from these equations and its evolution can be studied but the coupling 
between the equations makes such a method very difficult. Thus we adopt another approach, the one used in [10] for 
studying the symmetric Gaussian case. As shown in [9,10], the proper form of F-P equation can be achieved only 
after a further reparametrization, namely, t = — £!~ 2 lnCos(£!:r) and by studying the evolution in terms of r. This 
suggests us to use the same parametrization for our case also which gives 



with h = Vl - e" 2n2T and fh = e" a2r . 

Given an ensemble of the matrices V distributed with a probability density p(V) in the Hermitian matrix space, 
let P(p, t) be the probability of finding eigenvalues Xi[V] of H between m and pi + dpi at "time" r for an arbitrarily 
chosen H which can be expressed as follows, 



Note here that P(p,r) is in fact the conditional probaility P(p,r\p ) with p as the eigenvalue matrix of H but for 
simplification we use the former notation. As the r-dependence of P in eq.(2) enters only through Aj, a derivative of 
P with respect to t can be written as follows [10] 



II. GENERALIZED GAUSSIAN CASE 



H{t) = H fh + VhVL 
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The further treatment of above equation depends on the symmetry-classes of the matrices V and H, that is, whether 
they are real-symmetric or complex hermitian. 

In the real-symmetric case, can be expressed in terms of the matrix elements of the orthogonal matrix O which 
diagonalizes H, H = O t .A.O with A as the eigenvalue matrix (Appendix A). Using this in eq.(3) leads to the following 
form with gu = (1 + Ski), 
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an application of the partial integration on the second term in the right of above equation then gives 
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For simplification let us assume a same variance for all the diagonal matrix elements such that gkiau — ot for k = I. 



Now by expressing 
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and with the help of real-symmetric analog of 



relations (A3-A8) as well as the partial integration technique we can reduce eq.(5) as follows 
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By applying the orthogonality relation of matrix O, the eq.(6) can now be rewritten as follows 
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where (3=1 and 
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Following the similar steps and the appropriate relations given in Appendix A, one obtains the same equation in 
complex hermitian case too (with (3 = 2). The further analysis of the eq.(7) depends on the relative values of variances 
for the off-diagonal matrix elements. For a clear understanding, let us first consider a case where all the off-diagonals 
have the same variances, such that gkiotki = en' (for k ^ I). F then can be expressed in terms of the derivative of P 
with respect to y = 2L, 



F = 2h- 2 n 2 (l-y)y 



dP 
dy 



(9) 



Here we have used an idea quite important to our analysis that the variances of the matrix elements of 
V can also vary and therefore eigenvalues evolve not only as a function of the perturbation strength but the 
variances too. Note here that F can also be expressed as an extra drift term (F = J2 n d gjf" > ) , with R = 

I Tl,k<A9k^ a ki ~ aK7 dli ^(^i — K)O n kO n i) p(V)dV and it can be of interest to study the effect of this extra 
drift on the dynamics. For example, for R nonlinear in P, the motion may loose its random behaviour, showing some 
sort of periodicity. But in this paper, we restrict ourselves just to explore the universality of the dynamics, that is, to 
verify that P satisfies a similar F-P equation as the one proposed by Dyson. 

The eq.(7) contains the derivatives of P with respect to two parameters, namely, r and y. In order to prove that the 
statistical quantities in the GG case evolve in a similar way as in the SG case, one should be able to reduce the two 
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parameteric dependence of eq.(7) in a single parameter. We assume that it is possible along a curve parameterized 
by (j), t — t{4>), y = y{4>) such that 



d 



dr d dy d 

d<p dr dcf) dy 



(10) 



A comparison of eq.(10) with eq.(7) then leads us to the conditions for existence of such a curve namely, fr = 1 



and || = — y ^2 yS> ■ By solving these equations, one obtains the following curve in the (y, r)-space, 

e 2oV _ 1 



V 



e 2n 2 r _|_ y Q 



and r = cf> + r 



(11) 



where tq and j/o are arbitrary constants. Note the conditions mentioned above do not impose any constraints 
on these constants except that the curve should remain on the positive side of the upper-half of y — t plane. The 
extraction of tq from eq.(ll) leads to the following form of </> 
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where d = -^zn and can be chosen as unity (as the corresponding value of <f> still satisfies the required conditions). 
Thus, for y = 2L, eq.(ll) represents a set of curves along which eq.(7) adopts the form of a F-P equation 
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here the steady state is achieved for <fi — > oo which corresponds to t — > oo and a — > a'; the steady state solution is 
given by Yl iK j — ^ l^e - " S fc ^fc. Note that only r — > oo (with ^ = 0) no longer represents the steady state as 
in SG or WD case but represents a transition state with cf>(y) = — — l^log(rf -^-^ ). 

The steady state limit of P at (y,r) — (l,oo) seems to indicate that if the initial perturbation is such that the 
width (oc — ) of the eigenfunctions of H is smaller than the overlapping between them (oc ^7), they have a tendency 
to overcome this difference and extend over whole of the available Hilbcrt space. This also seems to suggest that a 
complicated interaction giving rise to larger transition probabiliies between energy levels as compared to the probability 
of staying in the same level can only be transient in nature. Given freedom to change, it rearranges itself in such a 
way so as to equalize these probabilities. 

This can be seen directly from eqs.(7,8) because in TV — > 00, £1 — > and H — > V thus removing any ^-dependence 
of A,'s and making F — > 0. A same result, for SG type perturbation in large iV-limit, has already been obtained [9]. 
The analogous behaviour in both GG and SG case can be attributed to the negligible effect of the distribution of 
A-diagonal matrix elements on the dynamics in presence of N(N — l)/2 off-diagonal ones. This indicates that, in the 
thermodynamic limit, difference of the variances does not leave its trace on the evolution of eigenvlaues. 

Again, in N — > 00 limit, all the three distributions, namely, SG, GG and WD evolve in the same way for arbitrary 
initial conditions. Thus all moments of the eigenvalues in SG and GG cases at any value of the transition parameter 
will be equal to the corresponding moments of particle positions of a WD gas undergoing Brownian motion due to 
the presence of thermal noise. However, as discussed in [9], this does not imply a Brownian motion of eigenvalues for 
the first two case. This is due to the difference of sources randomizing the motion of eigenvalues. For SG and GG 
cases, the source is the matrix V, acting like a quenched disorder while, for WD gas, the thermal noise gives rise to 
an annealed randomness. Although these different origins of randomness do not affect the static (for one parameter 
value) as well as the 2™ d -order parametric correlations but, as shown in ref.[9], the multiple parametric correlations 
higher than the 2 nd -order are different for the perturbed systems and WD gas. However, the origin of randomness 
being similar for the SG and GG cases, all the conclusions obtained for the correlations for the former [9] are also 
valid for the latter in the thermodynamic limit. 

It is interesting to note that eq.(4) can also be written in the following form (by first using the real-symmetric 
analog of eq.(A3) in eq.(4) and then taking the derivative with respect to /j, n inside the integral) 
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Now an application of the partial integration to the second term on the right of this equation gives 
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W klP (V)dV = JHS(fi t - Xi) [1 + 2a kl V 2 ] P (V)dV 



(15) 



By expressing this in terms of derivatives of P with respect to a, a' as before, we obtain a drift equation (without 
a diffusion term as well as the repulsive pairwise potential) 



dP(jH,r) 



d 
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where <j)\ describes a curve in a 3-dimensional (y, a, a') space such that = ^ 



2fT 
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da 
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above possibility of reduction of the F-P equation in a pure drift equation just by going to a higher parametric space 
implies the total blindness of eigenvalues to any mutually repulsive potential or locally fluctuating force in this space. 
Note that it is already well-known that a repulsion between the eigenvalues vanishes if the number of parameters of 
the hamiltonian undergoing variation is more than one in real-symmetric case (two for complex- hermitian case) [11]. 
However the same situation seems to exist even when the parameters undergoing variation are the ones governing the 
distribution of matrix elements of hamiltonians. 

So far we have considered the case where all the off-diagonal matrix elements Vij have same variance although 
different from that of diagonal. But in many physical situations e.g. in disordered systems, there occurs localization 
of eigenfunctions which may give rise to off-diagonal elements with different variances. One such situation most often 
encountered is where V^'s decay exponentially with respect to the distance from diagonal. It can be modeled by 
taking different variances for different off-diagonals, that is, guUki — ov, (r = \k — F m such a case is given 

by, with y r =^r 



F = 2h- 2 n 2 Y J ( 1 ~yr)vr^- (17) 

We assume that it is possible along a curve parameterized by <fi, r = t (</>), Y — Y(yi, .., yjv) such that 

d_ _ drd_ dY_d_ 



with -Jjy = J2 r 7>y 9y~ ■ B- easomn g again as before one finds that the universality of evolution of the probability 
density still survives but only by going to a higher (./V-dimensional) parametric space (r, y\, .., yjv), along a set of 
curves parametrized by <\> and given by the conditions = 1 ^ = 2h~ 2 il 2 with ^ = (1 — y r )y r - Again <p can be 
obtained by solving these equations, 



= T -2^ l0g 



Vr 



(e 2H-T _ 1} 



(19) 



Note here that J7 r contains the conribution only from those y r 's for which y r — 1 ^ 0. 

One can also consider a physical situation when all the off-diagonal matrix elements have different probability laws 
Quotki = a, (for k = I) and a k i ^ otij if (kl ^ ij). In this case, F now turns out to be as follows with y k i = 



dP 



F = 2h- 2 n 2 y2(l-y kl )y kl - 



(20) 



Now one has to go to N ^ N 2 ^ + 1-dimensional parametric space to recover the F-P equation similar to Wigner- 
Dyson gas. The curves in this case are again parameterized by (j> which is still related to r and Y by eq.(18) but now 
Y = Y({y k i}) is such that ^ = J2k<i ^5V"d§^ wu ^ n T$F = (1 ~ Vkl)ykl', the two other derivatives of t and Y with 
respect to <f) remain the same as in the above case. Proceeding exactly as before, <fi in this case can be shown to be 
following, 
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Again, as before, the Hfc<2 contains contributions from all yu 7^ 1. 

For SG case, <f> appearing in the correlators of the type < d(E, <j)) 1 d(E, 0) > involves variation of just one parameter, 
for CG case, two parameters or more. However, in both the cases, the expression of second order correlators is same 
due to same form of the F-P equation. Our study therefore suggests that the second order correlators involving 
variation of more than one parameter can always be expressed as those involving just one effective parameter. For 
example the correlators belonging to RM models of both localized as well as delocalized case can be expressed in the 
same form by using parameter (f> although the definition of <j> is different in two cases. 

Our study also suggests that nature of the dynamics of the eigenvalues is very sensitive to the number of parameters 
undergoing a change. For finite N, the motion may appear as pure drift suggesting no interaction between eigenvalues 
and a total absence of local fluctuating and frictional forces if variances of all the matrix elements are allowed to vary 
independently It seems Brownian as a function of just the relative variations of the variance. Under variation of only 
perturbation strength parameter, the presence of an exta drift like term makes the motion non-Brownian. But the 
contribution of this term being zero in N — > 00 limit, it recovers its Brownian nature. It can be of interest to study 
the effect of this extra drift on the dynamics. 

Although the study presented here deals with the energy level-dynamics of random matrix models of hcrmitian 
operators e.g. hamiltonian of complex systems, the results are also valid for the level-dynamics of unitary operators 
U e.g time-evolution operator being perturbed by a random matrix V taken from the GG ensembles, (U(t) = 
Uoexp[iy/rV] w Uo[l + iy/rV] for small r values). The required F-P equation in this case turns out to be the following 



dP(lH,(/>) _ 1 d 



d (3 \ 1 H m — fl n 

— h — > cot ' 



(22) 



with 4> defined as before. To prove the equivalence for U a symmetric unitary matrix ((3=1), one has to follow the 
similar steps as in the real-symmetric case of hermitian matrices as both are invariant under orthogonal transformation. 
Similarly for a general unitary U ((3 = 2), the steps are similar to the complex hermitian case. As obious from eq.(22), 
this is similar to the F-P equation obtained if V belonged to the SG ensembles [12]. . 



III. CALCULATION OF P(/x,t) 

In the preceeding section, we obtained a F-P equation governing the evolution of probability P([i, r) for an arbitrarily 
given H(<f>o) (the conditional probability). Here <po corresponds to the case r = and and H((f) ) can be obtained 
by using the relation H(<f>) = H((j)o)e~ n (*~^°) + V\/l — c~ 2 ° 2 (^~^°), substituting <j> in terms of r and y there and 
then comparing it with the expression for the hamiltonian H(t) given by eq.(l). Here V is a matrix with an invariant 
probability distribution (p(V) oc e~ aTrV ) and H(4>) can be written in the above form as it gives the correct evolution 
equation for the probability distribution. 

The formal similarity of the F-P equation with that of WD gas case as well as SG case makes it easier to obtain 
the P(fi, t) at least for (3 — 2 case as the solution for the latter case is already known. The P(fi, <fi) (and therefore 
P(H,t) can also be calculated directly from hamiltonian H(4>) by using sum of the matrices technique (that is, by 
evaluating a two- matrix integral) as for the SG case [4] because now both the component matrices have invariant 
distribution unlike eq.(l). For completeness purpose, We give here few of the steps, used in solving SG case, for our 
case; for details refer to [12]. 

As can readily be checked, both eq.(13) and (22) can also be written as follows (with a = 1 for simplification) 

where IQnI^ — |A(^)pe~ w Sf= Mfc with A (jit) = FIi<j (A 4 * ~ Mi) f° r the hermitian case and = ILk/c^ 11 { ^ 3 ^ k ^ for 
the unitary case. The transformation * = P/IQn] 13 ^ 2 allows us to cast eq.(23) in the suggestive form 

- = (24) 
where, for the hermitian case, the 'Hamiltonian' H turns out to be the Calogero-Moser Hamiltonian 



G 



Similarly, for the unitary case, H is the Calogero-Sutherland hamitonian 



H = E L + E cosec2 - H) - ^ N ( N " 1) (26) 

With the parabolic-confining potential (or periodic boundary conditions) and under the requirement (to take account 
of the singularity in H) that the solutions vanish as — Hj\^^ 2 when ^ and /x,- are close to each other, H, both in 
eq.(25) and (26), has well-defined (completely symmetric) eigenstates Cfe and eigenvalues Afc. This allows us to express 
the "state" VP and therefore P(n, <f>\Ho) as a sum over eigenvalues and eigenf unctions of H 



PMMfa)) = i^W^^^eM-W-MKkWCM (27) 



fc>0 



where /xo = (/ioi, M02, Mow) are the eigenvalues of H(<fio)- The joint probability distribution P(/x, (j>) can then be 
obtained by integrating over all initial conditions 

P(ji, 4>)= I P{H, 4>W MPW toWo (28) 



which further leads to correlations using standard techniques [12,13]. 

The eqs. (27-28) represent the formal solutions of eq.(13) (or eq.(22)). To proceed further, one needs to know the 
eigenvalues and eigenfunctions of H so as to express P in a compact form but these are explicitly known only for 
(3 — 2 case [8]. This is because for (3 = 2 the interaction term in eqs. (25-26) drops out and then P can be obtained 
explicitly. For CM model (eq.(25)), it is given as follows [4], 

P(fi,<j>\H(cf> );(3= 2) cx \^L\deb[f m (m - wy; <f> - <h)]ij=i...N (29) 



'AM' 



with f m (x -y;t) = cxp 



, and, for CS model [12], 
1 ■ A(/x) 



P(p, cj>\H(cj) y,p = 2) - -L|^S-|det[/ a (/ii - fJ-oj ; <j> - <h)]ij=i...N (30) 



where f s (x) = ± exp[-fc 2 + ikx]. 

As we already know P(/j,, <j>) and various correlations for WD case (for which <p = t) starting from various initial 
conditions [12,13], one can obtain these measures for GG case just by substituting <f) by its appropriate relationship 
with r and variances of the perturbation-matrix elements. For example, for the case H = Vil^ 1 (which corresponds 
to H(t — > 00) in eq.(l)) with different variances for different diagonals of V (eq.(17)), (j> = — (2fi 2 ) _1 log[]l^ =1 ^~^ ] 

(with d chosen to be unity for the same reason as in eq.(12)). Also note that, in this case, = 0. The initial 
condition <f> = now corresponds to all y r — > 00, thus implying the existence of only the diagonal elements of V 

and a Poisson distribution for P(fi ,4> = 0) (cx c ^ " with Vu = /xoi)- As (f> — > 00 when jy r — > 1 for all r, the 
equilibrium distribution is therefore given by the standard Gaussian ensembles. This case thus corresponds to the 
Poisson — > GE transition in the standard Gaussian ensembles [12,13] with </> now as a transition parameter (replacing 
t, the perturbation strength) with intermediate ensembles representing the cases for any choice of y r 's. The two-point 
correlations for this transition for V belonging to the SG ensemble (complex-hermitian) have already been obtained 
[13]. It should be noted that, as for the SG ensembles, <p must be rescaled to see the smooth transition [14]. . 



IV. POLYNOMIAL CASE 



In section II, we showed the equivalence of the dynamics of the Wigner-Dyson gas with the evolution of the 
eigenvalues of a hamiltonian under a perturbation V taken from a generalized Gaussian ensembles. However the 
claim about universality of dynamics and density correlators can only be made if the equivalence is proved, at 
least in thermodynamic limit, for a V with distribution of more general nature. In this section, We attempt to 
verify the claim by taking p(V) as the exponential of the polynomial form of V, p(V) = Ccxp(— ^2 k<[ Q(Vki)) 

with Q(x) = lki{r)x 2r (a polynomial of x with degree 2M), C as the normalization constant (referred as 
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polynomial case I later on) and variances for diagonal and off-diagonal matrix elements chosen to be arbitrary. Note 
the universality of the correlators for the case H = V, with V distributed as in this case, has already been studied in 
Ref. [14] (which corresponds to the steady state limit of the study given here) 

To obtain a F-P equation, we now need following equality which can be proved by using integration by parts, 



/[V] 



By using the real symmetric analog of eq.(Al) in eq.(3), followed by above equality, we get 
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p(V)dV 
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(32) 



Now for simplification let us consider a case where gkilki{r) = 7(r) if (& = and gkilkl{f) = l'{ r ) if 7^ 0- 
Using the tools given in appendix A and by an extensive use of the integration by parts while dealing with nasty 

integrals, the first term appearing in above equation can now be reduced in the same form as in GG case. This results 

in the following 



<9r 



7 '(1) ^ d^n 



d 



t^ + E 







P-z 



where (3—1 and 



n 2 1 



h 2 7(1) 7'(1) 
with F, G r i, G r 2 given as follows 



1 O 2 k 

/f1V h 2 ^ 



r=2 



7» r , ,l{r) n 
7 '(1) 7(1) 



^=E^/n^-^) 



Ea^(n^-Ai)^-s fcfc 



E ^^rt T/2r-l 



fc<Z 



dVkk 

p(v)dv 



p(y)dy 



E ^^rt y2r-l 
. dV kk kk 



p(V)dV 



We apply partial integration to F and rewrite it as follows, 

M 



p(V)AV + J2n(r)G r2 

r=2 



similarly G r \ and G r2 can be rewritten as follws, 



k<l J i 



M 



(2r - iK 2 r 2 - 2 7 '(i)nr - 2^ S7 '( S )^r 



+2s-2 



s=2 



p(V)dV 



(33) 



(34) 



(35) 
(36) 
(37) 



(38) 



(39) 



°- 2 = e/ n<^- a * 



(2r-l)F fc 2 r 2 



7(1)^ 



M 

E^(*K 2 r 2s - 2 

s=2 



p(V)dV 



(40) 
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For the complex Hermitian case also, one obtains an equation similar to eq.(33) with (3 = 2 and Z given by eq.(34). 
Note that all the terms appering in the expressions of F, G r \, GV2 are of the type J YliiHi — K)V J p(V)dV , and, as done 
in section II for GG case, they can easily be rewritten in terms of the derivatives with respect to 7(7-) 's if j < 2M. It 
is difficult to do so for terms with j > 2M, they probably could be expressed as higher order derivatives with respect 
to more than one variance parameters (this is the case at least for N = 2). However as physically significant RM 
models of complex systems generally correspond to limit TV — ► 00, it would be sufficient, for our purpose, to study the 
behaviour of the terms in this limit. As obvious from eqs. (38-40), the A or ^-dependence of F, G r \ and G r 2 is due to 
presence of the Aj's which can be clearly seen by using the equality S(X~ 1 [p\ — V) = 5(p — \[V])det\M^\ and rewriting 
these integral in terms of the function A _1 (^i). Let us consider the case when both 7(7*) and 7(7"'), r = 1 — > M, 
are independent of TV and all the matrix elements of V are distributed with zero mean; the latter choice implies the 
N-independence of Vki also. Thus the TV-dependence of F, G r \ and GV2 in this case is of the same order as appears in P 
and therefore it follows from eq.(34) that the contribution from Z to eq.(33) is negligible (as compared to the diffusion 
term and the drift term due to mutual repulsion) in that limit. Similarly for GG case, Z vanishes for A — > 00; Z here 
can be obtained by substituting 7(7") = 0, "f'(r) = 0, for r > 1, in cq.(13). Further for the case where distribution, 
although basis independent, contains a non-Gaussian term, that is, for p(V) oc exp[— a Tr(V 2 ) — X^2 7( r ) Tr{V 2r )\ 
(polynomial case II), the Z can be shown to be following, 

r=2 

p(V)dV^0 for TV ^00 (41) 

In thermodynamic limit (A — ► 00), the eq.(33) is same as the F-P equation governing the Brownian motion of 
particles in Wigner-Dyon gas (also the Sutherland model) with particle positions and time in the latter replaced 
by jUj and r in the former. This remains valid also for the case when Q(x) is chosen to be an arbitrary function 
expandable in a Taylor series. Thus we find that, under a perturbation taken from a ensemble with a sufficiently 
general distribution, the distribution of eigenvalues of the quantum system evolves in the same way as the distribution 
of particle position in the Wigner-Dyson gas for arbitrary initial conditions. This also implies that all the moments 
of the eigenvalues calculated for a parameter value r will be equal to the corresponding moments of particle positions 
of the latter. But the motion of the eigenvalues is not Brownian because the randomness in their motion comes from 
the matrix V which acts like quenched disorder [9]. 

This also implies the equivalence of second order parametric correlators in the two cases because they can be 
expressed as a sum over various moments of eigenvalues, weighted suitably and then averaged over equilibrium initial 
conditions. However, the F-P equation being Markovian in nature, its equivalence for the two cases can not lead us 
to a similar conclusion for the the higher order correlators which involve moments at more than one parameter value 
and therefore multiple parameter averaging. As mentioned in [9], the n-point correlations for the Wigner-Dyson gas 
can be expressed in terms of the two point functions, yielding an n-matrix integral while for RM models of quantum 
chaos, the n-point function remains a two matrix integral. 

It shuold be noted here that we have considered the case only for A-independent coefficients in the polynomial; a 
more careful analysis is required when these have different A-dependence. For example, if one or more coefficients 
are A-depcndcnt, increasing with TV increasing, the contribution to cq.(33) from terms in eqs. (38-40) is no longer 
negligible and therefore the evolution of P is no longer the same as in Gaussian case. 

V. CONCLUSION 

In this paper, we have analytically studied the response of energy levels of complex quantum systems to external 
perturbations modeled by generalized random matrix ensembles. Our results indicate the universality of two-point 
parametric density correlators as well as the static correlators of all orders thus agreeing with the numerically observed 
results for complex systems. One interesting feature revealed by our study is that for both the localized and the 
delocalized quantum dynamics of these systems, the second order parametric correlations can be shown to have the 
same form except that the definition of the effective parameter is different in the two cases. The method adopted 
here handicaps us from making any statement about universality of the higher order parametric correlations. 

Further, for the reasons given in the sections (I, IV), the correspondence shown between Wigner-Dyson gas and 
the quantum hamiltonian H = Hq + xV with V having a sufficiently general nature also implies the equivalence 
between the space-time correlators of the Calogero-Suthcrland system and the second order parametric correlations 



|f[<5(ui - Aj) 



E 



9kj 



d(v 2r - r ) 



kl 



dVki 



AaTr{V 2r ) - A^Tr{V 



2r+2s-2\ 
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of H. This further strengthens the idea contained in [1] and supported by studies in [9], namely, the ID Sutherland 
model provides a model hamiltonian for the dynamics of eigenvalues of quantum chaotic systems [10]. The equivalence 
shown between the F-P equation governing the evolution, of eigenvalues in GG and SG cases also turns out to be 
quite helpful in studying the correlations for various difficult but physically significant situations (depending on the 
variances of perturbation matrix) in the former case, by using the technique given in section III. Note that our result 
is quite general as we have shown the equivalence for arbitrary choice of variances. As already mentioned in section 
II, the results obtained in sections II, III and IV are also valid for the level-dynamics of unitary operators U. Note the 
analogy between the statistical properties of non-equilibrium Circular ensembles which model the eigenvalue spectra 
of unitary operators and non-equilibrium standard Gaussian ensembles has already been proved [12]. 

Our study still leaves many important questions unanswered. For example, universality or its absence among 
higher order correlations of complex systems and their similarities or differences with Wigner-Dyson gas is not fully 
understood. Further the results here have been obtained for explicit averaging over a generalized random perturbing 
potential. It would be interesting to know whether similar conclusions can also be obtained for the complex systems 
where ensemble averaging is not valid (e.g. Billiards) and eigenvalue statistics should be computed as an average over 
the N eigenvalues. The intution suggests a positive answer as the sufficiently complex systems are quite likely to be 
self-averaging. This implies that a single choice of matrices from a particular distribution are representative in the 
large TV-limit of the entire distribution and therefore an average over eigenvalues should give the same result as an 
ensemble average. 
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APPENDIX A: 

For H a complex hermitian matrix, its eigenvalue equation is given by A = U.H.U + with A as eigenvalue matrix 
and U the eigenvector matrix which is unitary. Now for a complex hermitian perturbation Vij = Vy ; i + iViy.i = V* { 
(with Vu-2 = 0), the rate of change of eigenvalues and eigenvectors, with respect to various components of matrix 
elements of V, can be described as follows: 

= - fi2A « + f E ^7 [VirAUmU* nj + U nj U* ni ) + iV^U^ U nj U* m )] (Al) 

i<j y%j 



where 



Further 



oh ~ T vn , A „. 

— = -SPH + — (A2) 



" X " ■■ l! ' r) p nk K l + U nl U* nk } (A3) 



dV k i-i Q,g k i 
d\ n 1i(t) 



dVki-2 Qgu 



[U nk U* nl - U nl U* k ] (k I) (A4) 



= (if k = l) (A5) 



and 



StT = ^7 E T-^-U mp (U^ k U nl + U M (A6) 



DVki-i Qgki ^—f A m - A n 
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dV k i a 



-Mr) 



E 



U mp(U^ k U n l - U* ml U nk ) 



and 



E 

k<l 



OVki-i ' dVki-2 



2h(r) 



E 



(A7) 



(A8) 



For the real-symmetric case, the corresponding relations can be obtained by using Uij = U*j (as eigenvector matrix 
is now real-orthogonal) in cqs.(Al-A8) and taking V%y,i = for all values of i,j (can also be found in Ref.[9]). 



[1] B.D.Simons and B.L.Altshuler, Phys. Rev. Lett. 70, 4063 (1993). 

B.D.Simons, P.A.Lee and B.L.Altshuler, Phys. Rev. Lett. 70, 4122 (1993). 
[2] C. Beenakker, Phys. rev. Lett. 70 (1993). 
[3] K.B.Efetov, Adv. in Phys. 32, 53 (1983). 

[4] M.L.Mehta Random Matrice, (Acedemic, New York, 1991). 
[5] M.V.Berry, Proc. R. Soc A, 400, (1984). 
[6] F.Dyson, J. Math. Phys. 3, 1191 (1962). 

[7] M.Faas, B.D.Simons, X.Zotos and B.L.Altshuler, Phys. Rev. B 48, 5439 (1993). 

[8] B.Sutherland, J. Math. Phys. 12, 246 (1971); 12, 251, (1971); Phys. Rev. A 4 2019 (1971); 5, 1372 (1972). 
[9] O. Narayan and B.S.Shastry, Phys. Rev. Lett. 71, 2106 (1993). 
[10] B.S.Shastry, Proceedings of the 16th Taniguchi International Symposium on the Theory of 

Condensed Matter: "Correlation Effects in Low Dimensional Electron Systems", 

October 23 to 29, 1993, Shima, Japan (Ed. N.Kawakami and A.Okiji, Springer Verlag, 1994). 
[11] F.Haake, Quantum Signatures of Chaos, (Springer, Berlin, 1991). 
[12] A.Pandey and P.Shukla, J. Phys. A, (1991). 
[13] A.Pandey, Chaos, Soliton and Fractals, 5, (1995). 
[14] P.Shukla, Submitted to Phys. rev. Lett. (1998). 
[15] E.Brezin and A. Zee, Nucl. Phys. B402 613 (1993). 
[16] G.Hackenbroich and H.A.Weidenmuller, Phys. Rev. lett. 74 (1995). 
[17] Present Address: Department of Physics, 

Indian Institute of Technology, 

Kharagpur, West Bengal, India 

E-Mail : Shukla@phy.iitkgp.emet. in 



11 



